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Abstract 

Using a generalization of vector calculus for space with non-integer dimension, we con¬ 
sider elastic properties of fractal materials. Fractal materials are described by continuum 
models with non-integer dimensional space. A generalization of elasticity equations for 
non-integer dimensional space, and its solutions for equilibrium case of fractal materials 
are suggested. Elasticity problems for fractal hollow ball and cylindrical fractal elastic 
pipe with inside and outside pressures, for rotating cylindrical fractal pipe, for gradient 
elasticity and thermoelasticity of fractal materials are solved. 

PACS: 62.20.Dc; 81.40.Jj; 45.10.Hj 


1 Introduction 

Fractals are measurable metric sets with non-integer Hausdorff dimensions DEI' The main 
characteristic of fractal set is non-integer Hausdorff dimension that should be observed on all 
scales. The Hausdorff dimension is a local property, i.e. this dimension characterize (measure) 
property of a set of distributed points in the limit of a vanishing diameter, which is used to 
cover subset of the points. By definition the Hausdorff dimension requires the diameter of the 
covering sets to vanish. In general, real materials have a characteristic smallest length scale 
Ro such as the radius of a particle such as atom or molecule. In fractal materials the fractal 
structure cannot be observed on all scales but only those for which R > R 0 , where R 0 is the 
characteristic scale of the particles. For real materials, a non-integer mass dimension can be 
used instead of Hausdorff dimension. The mass dimension described how the mass of a medium 
region scales with the size of this region, where we assume unchanged density. For many cases, 
we have an asymptotic relation between the mass M(W) of a ball region W of material, and the 
radius R of this ball. The mass of fractal material satisfies a power-law relation M(W) ~ R D . 
The parameter D is called the non-integer mass dimension of fractal material. This parameter 
does not depend on the shape of the region W, or on whether the packing of sphere of radius 
Ro is close packing, a random packing or a porous packing with a uniform distribution of holes. 
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Therefore fractal material can be considered as a medium with non-integer mass dimension. 
Although, the non-integer dimension does not reflect completely the geometrical and dynamical 
properties of the fractal materials, it nevertheless permits a number of important conclusions 
about the behavior of materials. It allows us to use effective models that take into account 
non-integer dimensions. 

We can distinguish the following approaches to formulate models of fractal materials: 

1) Approach based on methods of ” Analysis on fractals” [3], 0| El E, 7, E] can be considered 
as the most rigorous approach to describe fractal materials. Unfortunately a possibility of 
application of the ’’Analysis on fractals” to solve real problems of fractal material now is very 
limited due to weak development of this area of mathematics. 

2) To describe fractal material we can apply special continuum models suggested in 0 
nm nu ng and then developed in the works [ 13] - [2D] and [21]. These models can be called 
the fractional-integral continuum models. In this approach we use integrations of non-integer 
orders, and two different notions such as density of states and distribution function [31]. The 
order of the fractional integrals is equal to mass dimension of fractal materials. The kernels of 
these integrals are defined by the power-law type of density of states. 

3) Fractional derivatives of non-integer orders are used to describe some properties of fractal 
materials. This approach has been suggested in papers [33] 23] [23], where so-called local 
fractional derivatives are used, and then developed in the works [25] [26] [37, 28] [29, [30] 131]. 
These models can be called the fractional-differential models. 

4) Fractal materials can be described by using the theory of integration and differentiation 
for non-integer dimensional spaces [32] [33] 33] • Fractal materials are described as continuum in 
non-integer dimensional spaces. The dimension of the spaces are equal to the mass dimensions 
of fractal materials. 

Unfortunately there are not enough differential equations that are solved for various prob¬ 
lems for fractal materials in the framework of the fractional-differential model and by methods 
of ’’Analysis on fractals”. 

The fractional-integral continuum models are used to solve differential equations for various 
problems of elasticity of fractal materials B31031021 [El EOT, and thermoelasticity of fractal 
materials [70] E]. 

Continuum models with non-integer dimensional spaces are not currently used to describe 
elasticity of fractal materials. In this paper, we consider approach based on the non-integer 
dimensional space to describe elasticity of isotropic fractal materials. The main difference of 
the continuum models with non-integer dimensional spaces and fractional-integral continuum 
models suggested in 0 m ei eh may be reduced to the following: (a) Arbitrariness 
in the choice of the numerical factor in the density of states is fixed by the equation of the 
volume of non-integer dimensional ball region, (b) In the fractional-integral continuum models 
the differentiations are integer orders whereas the integrations are non-integer orders. In the 
continuum models with non-integer dimensional spaces the integrations and differentiations are 
defined for the spaces with non-integer dimensions. The power law M ~ R D can be naturally 
derived by using the integrations in non-integer dimensional space [33], where the dimension of 


2 


this space is equal to the mass dimension of fractal material. 

A vector calculus for non-integer dimensional space proposed in this paper allows us to 
use continuum models, which are based on non-integer dimensional space, to describe fractal 
materials. This is due to the fact that although the non-integer dimension does not reflect all 
geometrical and dynamical properties of the fractal materials, it nevertheless allows us to get 
important results about the behavior of fractal materials. Therefore continuum models with 
non-integer dimensional spaces can describe a wide class of fractal materials. 

Integration over non-integer dimensional spaces are actively used in the theory of critical 
phenomena and phase transitions in statistical physics [35, 36], and in the dimensional reg¬ 
ularization of ultraviolet divergences in quantum field theory [33 |38l [32]. The axioms for 
integrations in non-integer dimensional space are proposed in [3D] [33] and this type of inte¬ 
gration is considered in the book by Collins [32] for rotationally covariant functions. In the 
paper [33] a mathematical basis of integration on non-integer dimensional space is given, and 
a generalization of the Laplace operator for non-integer dimensional spaces is suggested. Using 
a product measure approach, the Stillinger’s methods [33] has been generalized by Palmer and 
Stavrinou [33] for multiple variables case with different degrees of confinement in orthogonal di¬ 
rections. The scalar Laplace operators suggested by Stillinger in [33] and by Palmer, Stavrinou 
in [31] for non-integer dimensional spaces, have successfully been used for effective descriptions 
in physics and mechanics. The Stillinger’s form of Laplacian for the Schrodinger equation in 
non-integer dimensional space is used by He [ID] E, 42] to describe a measure of the anisotropy 
and confinement by the effective non-integer dimensions. Quantum mechanical models with 
non-integer (fractional) dimensional space has been discussed in [33] 33] 03, 3H 05] 05] 37] 
and mi mmm- Recent progress in non-integer dimensional space approach also includes 
description of the fractional diffusion processes in non-integer dimensional space in [52], and 
the electromagnetic fields in non-integer dimensional space in [53] [53, 55] and [55! [571 58, [59]. 

Unfortunately, in the articles [33] [33] are proposed only the second order differential op¬ 
erators for scalar fields in the form of the scalar Laplacian for the non-integer dimensional 
space. A generalization of the vector Laplacian [5T] for the non-integer dimensional space is 
not suggested. The first order operators such as gradient, divergence, curl operators, and the 
vector Laplacian are not considered in [33] [33] also. In the work [60] the gradient, divergence, 
and curl operators are suggested only as approximations of the square of the Laplace opera¬ 
tor. Consideration only the scalar Laplacian in non-integer dimensional space approach greatly 
restricts us in application of continuum models with non-integer dimensional space for fractal 
materials and material. For example, we cannot use the Stillinger’s form of Laplacian for the 
displacement vector field u(r, t) in theory of elasticity and thermoelasticity of fractal materials, 
for the velocity vector field v(r, t) in hydrodynamics of fractal fluids, for electric and magnetic 
vector fields in electrodynamics of fractal media in the framework non-integer dimensional space 
approach. 

In this paper, we define the first and second orders differential vector operations such as 
gradient, divergence, the scalar and vector Laplace operators for non-integer dimensional space. 
In order to derive the vector differential operators in non-integer dimensional space we use 
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the method of analytic continuation in dimension. For simplification we consider rotationally 
covariant scalar and vector functions that are independent of angles. It allows us to reduce 
differential equations in non-integer dimensional space to ordinary differential equations with 
respect to r. The proposed operators allows us to describe fractal materials to describe processes 
in the framework of continuum models with non-integer dimensional spaces. In this paper we 
solve elasticity problems for fractal hollow ball with inside and outside pressures, for cylindrical 
fractal elastic pipe with inside and outside pressures, and for rotating cylindrical fractal pipe, 
for gradient elasticity and thermoelasticity of fractal materials. 

2 Differential and integral operators in non-integer di¬ 
mensional space 

To derive equations for vector differential operators in non-integer dimensional space, we use 
equations for the differential operators in the spherical (and cylindrical) coordinates in M n for 
arbitrary n to highlight the explicit relations with dimension n. Then the vector differential 
operators for non-integer dimension D can be defined by continuation in dimension from integer 
n to non-integer D. To simplify we will consider only scalar fields tp and vector fields u that 
are independent of angles 


<p(r) = u(r) = u(r) = u r e r , 

where r = |r| is the radial distance, e r = r/r is the local orthogonal unit vector in the directions 
of increasing r, and u r = u r (r ) is the radial component of u. We will work with rotationally 
covariant functions only. This simplification is analogous to the simplification for definition of 
integration over non-integer dimensional space described in Section 4 of the book [32]. 


2.1 Vector differential operators for spherical and cylindrical cases 

Explicit definitions of differential operators for non-integer dimensional space can be obtained 
by using continuation from integer n to arbitrary non-integer D. We mote that the same 
expressions can be obtained by using the integration in non-integer dimensional space and the 
correspondent Gauss’s theorem. 

We define the differential vector operations such as gradient, divergence, the scalar and 
vector Laplacian for non-integer dimensional space. For simplifications, we assume that the 
vector field u = u(r) be radially directed and the scalar and vector fields </?(r), u(r) are not 
dependent on the angles. 

The divergence in non-integer dimensional space for the vector held u = u(r) is 


Divf u 



D - 1 
- u r . 

r 


( 1 ) 
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The gradient in non-integer dimensional space for the scalar field p = p{r) is 


Gradf p — ^ e r . 
or 


( 2 ) 


The scalar Laplacian in non-integer dimensional space for the scalar field p = p(r ) is 


S aD,„ _ /"i- jD , D — 1 dp 


’ A"^ = Div" Grad" p = 


+ 


The vector Laplacian in non-integer dimensional space for the vector field u = M r (r) e r is 


(3) 


Va-D 


A"u = Gradf Divf u = ^ 


<9 2 w r D — 1 chq 


+ 


D 


r dr 


u r e r . 


(4) 


If D = n, equations ([IMl) give the well-known formulas for integer dimensional space M n . 
We can consider materials with axial symmetry, where the fields p{r) and u(r) = u r (r)e r 
are also axially symmetric. Let the Z-axis be directed along the axis of symmetry. Therefore 
we use a cylindrical coordinate system. 

The divergence in non-integer dimensional space for the vector field u = u(r) is 


D du r D - 2 
Div" u = H- u r . 


dr r 

The gradient in non-integer dimensional space for the scalar field p = p{r) is 


(5) 


Gradf p = 


dp 

dr 


( 6 ) 


The scalar Laplacian in non-integer dimensional space for the scalar field p = p{r) is 

s a d _ d 2 p D- 2 dp 
r * dr 2 r dr ' 

The vector Laplacian in non-integer dimensional space for the vector field u = v(r) e r is 


(7) 


VaD 


A"u = 


d 2 u r D — 2 du r D — 2 


dr 2 


dr 


■ u r ] e 


( 8 ) 


Equations ([MS]) can be easy generalized for the case p = p(r, z ) and u(r, z) = u r (r , z ) e r + 
u r (r, z) e z . In this case the curl operator for u(r, z) is different from zero, and 


Curl" u = 


du r dv. 


dz dr 


e*. 


(9) 


Equations (JTj ) - (J8| ) with D — 3 and (JUJ) give the well-known expressions for the gradient, 
divergence, curl operator, scalar and vector Laplacian operators 
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The proposed operators for 0 < D < 3 allows us to reduce non-integer dimensional vector 
differentiations (JTJ) - ([4]) and (J5]) - (JEJ) to derivatives with respect to r = |r|. It allows us to reduce 
partial differential equations for fields in non-integer dimensional space to ordinary differential 
equations with respect to r. 

For a function <p = ip(r, 6) of radial distance r and related angle 6 measured relative to 
an axis passing through the origin, the scalar Laplacian in a non-integer dimensional space 
proposed by Stillinger [33] is 


St^D 


1 d ( D _ 1 d \ 1 d ( . D _ 2 d \ 

-'-dr h ~ Tr) + r 2 '*36) ’ 


( 10 ) 


where D is the dimension of space (0 < D < 3), and the variables r > 0, 0 < 6 < n. Note that 
( s *A d ) 2 ^ st A 2D . If the function depends on radial distance r only [ip = then 


St A D ip(r) 


1 9 Ld -i \ 

r D ~i Q r dr J 


d 2 ip{r ) D « 1 dip{r) 


( 11 ) 


It is easy to see that the Stillinger’s form of Laplacian St A D for radial scalar functions ip(r) = 
cp(r) coincides with the scalar Laplacian s Af defined by (J3J), i.e., 

St A D ip(r) = S A D p>{r). (12) 


The Stillinger’s Laplacian can be applied for scalar fields only. It cannot be used to describe 
vector fields u = u r (r) e r because this Laplacian for D = 3 is not equal to the usual vector 
Laplacian for M 3 , 

st A 3 u(r) ^ Au(r) = - 4 Ur ) ( 13 ) 

' OT ^ T C/l™ y 1 " j 

The gradient, divergence, curl operator and vector Laplacian are not considered by Stillinger 
in paper [33] . 


2.2 Integration over non-integer dimensional space 

Integration for non-integer values of dimension D is defined by continuation in D [38,132]. The 
following properties suggested in |39j for integrals in D -dimensional space are necessary for 
applications [32] : 

a) Linearity: 

J (a/i(r) + 6/2 ( r)) d D r = a J /i(r) d D r + b J / 2 ( r) d D r, (14) 

where a and b are arbitrary real numbers. 

b) Translational invariance: 

J f(r + r 0 )d D r= j f(r)d D r (15) 
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for any vector r 0 . 

c) Scaling property: 

J f{Xr)d D r = X~ D J f(r)d D r (16) 

for any positive A. 

Note that linearity is true of any integration, while translation and rotation invariance are 
basic properties of a Euclidean space. The scaling property embodies the ^-dimensionality. 
Not only the above three axioms are necessary, but they also ensure that integration is unique, 
aside from an overall normalization [39], These properties must be used in order to have non¬ 
integer dimensional integrations [32]. These properties are natural in application of dimensional 
regularization to quantum held theory [32] . 

In general, we can consider any functions of the components of its vector argument r. 
However, we do not know the meaning of the components of a vector in non-integer dimensions. 
In this paper, we will work with rotationally covariant functions for simplification. So we will 
assume that / is a scalar or vector function only of scalar products of vectors or of length of 
vectors. For example, in the elasticity theory, we consider the case, where the displacement 
vector u(r), is independent of the angles u(r) = u(r), where r = |r|. The integration defined 
by equation (TT7l) satisfies the properties (TT4|) - (fT6lh 

The non-integer dimensional integration for scalar functions /(r) = /(|r|) can be defined in 
terms of ordinary integration by the equation 


/ 



dVlD—i 


dr r 


D—l 


f(r), 


(17) 


where 



2tt d / 2 

T(D/2) 


Sd- i. 


(18) 


Equation (1T51) with integer D — n gives the well-known area SA-i of (n — l)-sphere with unit 
radius. 

As a result, we have [32] the explicit definition of the continuation of integration from integer 
n to arbitrary fractional D in the form 


J d ° r /(M) 


27^/2 

rp/2) 


dr 


„D-1 


/(f) 


(19) 


This equation reduced non-integer dimensional integration to ordinary integration. Therefore 
the linearity and translation invariance follow from linearity and translation invariance of ordi¬ 
nary integration. The scaling and rotation covariance are explicit properties of the definition. 

As an example of applications of equation (fl9|i . we can consider non-integer dimensional 
integration for the function 


/(r 2 ) 


r 2 + a 
r 2 + 6’ 


( 20 ) 
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where a and b are real numbers. The integral for (T20l) can be explicitly computed 


/ d ° r = ^ D ' 2 (“/ 6 - 1 ) r (i - o / 2 ). 

The other example is the non-integer dimensional integration is 

[ ri D r r2a _ + D/2) T(/3 — a — D/2) D / 2 D+2a _ 2 p 

J (r 2 + a?Y — T{D/2)T{0) 


(21) 


( 22 ) 


where r = |r|. 

3 Mass and moment of inertia for fractal materials 

3.1 Mass of fractal materials 

Fractal materials can be characterized by the relation between the mass M(W ) of a region W 
of fractal material, and the size R of the region containing this mass: 

/ TD \ d 

M d (W) = M 0 (- , Rj R 0 » 1. (23) 

The parameter D is called the non-integer mass dimension of fractal material. The parameter 
D, does not depend on the shape of the region W, or on whether the packing of sphere of radius 
Rq is close packing, a random packing or a porous packing with a uniform distribution of holes. 
The cornerstone of fractal materials is the non-integer mass dimension. The mass dimension of 
real fractal materials can be measured by box-counting method, which means drawing a box of 
size R and counting the mass inside. 

The fractality of material means than the mass of the region W C M 3 increases more slowly 
than the 3-dimensional volume of this region. For the ball region of the fractal medium, this 
property can be described by the power law Mr>{W) ~ R D , where R is the radius of the ball. 

Fractal material is called homogeneous if the power law M D {W) ~ R D does not depend on 
the translation of the region. The homogeneity property of the material can be formulated in 
the form: For all two regions W\ and W 2 of the homogeneous fractal material with the equal 
volumes Vo(fFi) = Vd(W 2 ), the masses of these regions are equal M^iWi) = Md(W 2 ). 

The power law (f23i) can be naturally derived by using the integration in the non-integer 
dimensional space such that the space dimensions is equal to the mass dimension of the material. 

The mass of the region W of fractal material in W can be calculated by the integral in 
non-integer dimensional space 

= [ P(r)d D r, 

Jw 
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M d (W) 


(24) 





where r is dimensionless vector variable. For a ball region W with radius R and density 
p( r) = po = const, we get the mass is defined by 


Mo{W) = Po Vo = r( ^ R°. (25) 

This equation define the mass of the fractal homogeneous ball. For D = 3, equation (125]) gives 
the well-known equation for mass of non-fractal ball M 3 = (Ap^) / ?>R i because T(3/2) = \/n/2 
and T(z + 1 ) = z T(z). 

3.2 Moment of inertia of fractal materials 

Let us consider a calculation of scalar moment of inertia I(t), which is used when the axis 
of rotation is known. The scalar moment of inertia of a rigid body with density p'(r',t) with 
respect to a given axis is defined by the volume integral 

m= [ pV, t) r'jL 2 dV, (26) 

Jw 

where (r')j_ is the perpendicular distance from the axis of rotation, and dV^ = dx' x dx' 2 dx' z . We 
note that SI units of I'ki is kg ■ m 2 . 

To generalize equation for non-integer dimensional space, we should represent this 
equation through the dimensionless coordinate variables. We can introduce the dimensionless 
values Xk = x' k / Ro, r = r'/R 0 , where Rq is a characteristic scale, and the density p(r,t) = 
Rq p'(r R 0 , t). SI units of p is kg. We define the following moments of inertia I(t) = RQ 2 I'(t). 
As a result, we obtain 

1= [ p(r) ri d 3 r. (27) 

Jw 3 

where x k (k = 1, 2, 3) and r are dimensionless. We note that SI units of Ri is kg. 

This representation allows us to generalize equation of the scalar moment of inertia for a 
fractal material 

I (D \t) = [ p(r,t) v\ dV D , (28) 

Jw 

where I) is a mass dimension of fractal material. 

3.3 Moment of inertia of fractal solid ball 

Let us consider a fractal solid ball with radius R , and mass M. Note that the component of 
the radius perpendicular is 

r l = ( r sin#) 2 . 
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Using the integration in a non-integer dimensional space, we have 


I (D) = / dr (r sin 9) 2 p(r, 6) = 


27t ( d_1)/2 r°° 


Jw r((L> - l)/ 2 ) J 0 

For homogeneous materials p(r) = p 0 , we get 


drr 


D -1 


d6 sin u 2 6 (r sin #) 2 p(r, 0). 


(29) 


/du = 


2vr( D - 1 )/ 2 


Po 


r ((-D - l)/ 2 ) Jo 


dr r' D+1 / d6 sin D 0 = 


2 tt ^- 1 )/ 2 po vt 1 / 2 r ( D /2) i ? D+2 


27T D / 2 r(L)/2)po i? D+2 


r((D - l)/2) r(£>/2 + 1) D + 2 

vr ^ 2 (D-l)po 


r(D /2 - 1/2) r(£>/2 + 1) D + 2 (£> + 2) T(D/2 + 1 ) 


i? D+2 , 


where we use 


d6 sin D 0 = 


7r 1 / 2 r(D/2) 

r(D/2 + i)' 


Using the expression for mass f[25lh we can rewrite 


as 


/(») = 


2r(D/2) 


(£> + 2) r(D/2 - 1/2) 


M d B 2 = ——1 Mr, R 2 , 
D + 2 


( 30 ) 

( 31 ) 

( 32 ) 


where we use T(z) = (z — l)r(z — 1 ). For D — 3 , equation (j 32 j) gives the well-known equation 
for moment of inertia of non-fractal ball / 3 = (2/5 )MK 2 . 


4 Elasticity theory of fractal material 

4.1 Elasticity theory of non-fractal material 

The linear elastic constitutive relations for isotropic case is the well-known Hooke’s law that 
has the form 

cr, 3 — Ac kkdij -f- “2fiEjj , (55) 

where A and p are the Lame coefficients, o VJ is the stress, £ki is the strain tensor. This expression 
determines the stress tensor in terms of the strain tensor for an isotropic material. 

For homogenous and isotropic materials, the constitutive relation (133|) gives the equation 
for the displacement vector fields u = u(r, t) in the form 

A graddivu + 2p A u + f = pD^u, (34) 

where f = u(r, t) is the vector field of external force density. 

If the deformation in the material is described by u(r, t) = u(r, t) e r , then equation ()34l) has 
the form 

(A + 2 n) Au(r,f) + f(r, t) = pD 2 u(r,t). (35) 
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A formal generalization of equations fl35l) for fractal material in the framework of non-integer 
dimensional models is 


(A + 2 ji) V A® u(r, t) + f (r, t) = p D^u(r, t ), (36) 

where v A^ is defined by (j4]) . Equation (1361) described dynamics of displacement vector for 
fractal clastic materials. 


4.2 Strain and stress in non-integer dimensional space 

Any deformation can be represented as the sum of a pure shear and a hydrostatic compression. 
To do so for fractal materials, we can use the identity 

&kl (&kl T “yrAfc/Sjj. (37) 


The first term on the right is a pure shear, since the sum of diagonal terms is zero. Here we use 
the equation 8a = D for non-integer dimensional space (for details see Property 4 in Section 
4.3 of [32]). The second term is a hydrostatic compression. For D = 3, equation (1371) has the 
well-known form 

l x 

+ —Okl £iii 

where 8 tl = 3 is used. 

The stress tensor can be represented as 

&kl A 8kl T 2/1 ^'£kl ~J^8f c i EaJ , (39) 



£kl ~ £kl 


r 8ki u 


where the bulk modulus (modulus of hydrostatic compression) is related to the Lame coefficients 
by 

K = A + A (40) 

In the hydrostatic compression of a body, the stress tensor is 

<7ki = -phi■ (41) 

Hence we have 

CT kk = -pD. (42) 

If we use the Hooke’s law for isotropic case in the form, then 

o'u — (A D + 2 fi) £jj. (43) 

The components of the strain tensor is 

I 

£ rr = -rr~ = (e n Grad ® u r ). (44) 

or 
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Using ([I|), and the trace of the strain tensor 

e(r ) = Tr[e k i} = £ kk = Div" u = — H--— u r . (45) 

Note that we can consider 

e(r) — e rr (r) = Divf 3 u — (e r , Gradf 5 u r ) =- u r 

r 

as a sum of the angle diagonal components in the spherical coordinates of the strain tensor. 
For D = 3, we have the well-known sum of the angle diagonal components in the spherical 
coordinates of the strain tensor 

_ 2 

EGO "F Emm Zl r . 

r 

When we consider the fractal medium distributed in three-dimensional space we can define the 
effective value of the diagonal angular components of the strain tensor 


F e ff — F e ff — 
L ee ~ t-cptp ~ 


D — 1 


2 r 


■ u r 


(46) 


Using (1441) and 
coordinates is 


, the components of the stress tensor a k i = cr k i(r) in the spherical 


du r D — 1 

a rr (r) = 2 n£ rr (r) + X e(r) = (2 /i + A) —-b A- u r 


(47) 


dr r 

ft is well-known that the diagonal angular components of the stress tensor for D = 3 in 
spherical coordinates are 


a ee {r) = 2/i£ eg (r) + Ae(r), cr w (r) = 2/ze w (r) + Ae(r). (48) 

For the fractal medium distributed in three-dimensional space we can define the effective 
value of the diagonal angular components of the stress tensor 

<?eo f ( r ) = 2 ^ £ ee f ( r ) + Ae(r), (49) 

rfj (r) = 2 [i £$(r) + A e(r). (50) 

Using the form of the effective components (1461) . we get 

a,, n_i 

= </(r) = A ^ + (A + f,) —— u„ (51) 


This equation define the diagonal angular components of the stress tensor for spherical coordi¬ 
nates. 
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4.3 Equilibrium equation for fractal materials 

For static case, equation (1361) . we have 

V A? u(r) + (A + 2p) _1 f(r) = 0, (52) 


where u = u r e r and f = f(r)e r . Here A and p are the Lame coefficients. Using (Hj), we 
represent equation (1521) in the form 


d 2 u r (r) D — 1 du r (r) 


D - 1 


u r {r ) + (A + 2/i) 1 f{r) = 0. 


The solution of (fool) is 

u r {r ) = Ci r + C 2 r 1_D + I f {D,r), 

where Ci and C 2 are constants dehned by boundary conditions, and 


(53) 

(54) 


If(D, r) = 

For f(r) = / 0 , we get 


D (A + 2/i) 
I f (D,r) = 


and the displacement is 


u r (r) = Ci r + C 2 r 1 u — 


J drr D f(r) — J drrf(r)J. 

_ for 2 

2 (-D + 1) (A + 2/i) 

for 2 


(D + 1) (A + 2/i) 


(55) 


(56) 


(57) 


The components of the stress tensor a ki = a k fr) for the spherical coordinates can be calculated 
by equations (UTD and (HTTP . 


4.4 Elasticity of fractal hollow ball with inside and outside pressures 

Let us determine the deformation of a hollow fractal ball with internal radius R\ and external 
radius R 2 with the pressure pi inside and the pressure p 2 outside. 

We can use the spherical polar coordinates with the origin at the center of the ball. The 
displacement vector u is everywhere radial, and it is a function of r = |r| alone. Then the 
equilibrium equation for fractal ball is 

(A + 2/i) u(r) = 0, (58) 


where u = u r e r . Using (J4J) , we represent equation (l52|) in the form 


d 2 u r (r) D — 1 du r (r) 

dr 2 r dr 


D — 1 


■ u. 


(r) = 0. 


(59) 
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(60) 


The solution of (1531) is 


u r {r) = Ci r + C 2 r 1 D 


The constants C\ and C 2 are determined from the boundary conditions for radial stress 


<J rr (yR\ ) Pi) (J rr (R 2 ) P2- 

Using that the radial components of the stress is 

du r D — 1 

CTrr(r) = (2 H + A) —-h A - U r 

or r 


we get 

Then we have 


a rr (r ) = (2 n + D A) Ci + 2 (1 — T>) fi C 2 r 


-D 


(2 p + D A) Ci + 2(1 — D) p R^ D C 2 — —Pi, 

(2 pL + D A) Ci + 2(1 — D) p D C 2 = —P2- 
As a result, the coefficients have the form 

-(p 1 R^ D ~P2Ri D ) -(P 2 R 2 -P 1 R 1 ) 


C! = 


Co = 


(2p + DX) (R^ d - R^ d ) (2/j + DA) (A 2 d - i2f)’ 

-(ps-pO^! A 2 ) d 


P2 “Pi 


2 (1 - D) n {R^ d - Ri D ) 2(1 -D)(i(R?-R?y 

Then the radial components of the stress is 


_ ^ -faRz-PiR?) , ( P 2-Pi)(RiR2r 

a rr\ r ) nfl nO ' l_ r>D tdD r 

1L2 Xti it 2 1 X 1 

The stress distribution in a ball with pressure pi—p inside and p 2 = 0 outside is gives by 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 


Orr r = 


pR[ 


D 


V 


■r~ D = 


p R\ 


D 


jdD tdD dD pD p D p D 

1L2 jILi lt2 it2 ./Xi 


1-1 + 




(69) 


The stress distribution in an infinite elastic material with spherical cavity with radius R 
subjected to hydrostatic compression 


a rr {r) = -p 



( 70 ) 


that can be obtained by putting R 1 = R, R 2 —* 00 , p\ = 0 and p 2 = p in equation (153|) . 
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5 Elasticity of fractal material with radial distribution 
in cylinder and pipe 

5.1 Elasticity of fractal pipe and cylinder 

If we use the cylindrical coordinates with the z-axis and the vector field u(r, t) is a purely radial 

u = u r (r) e r , (71) 


where e r = r/r, then 




d 2 u r {r ) D — 2 du r 
dr 2 r dr 


D -2 



6 r . 


(72) 


Note that we have (D — 2) instead of (D — 1) in this equation. For D — 3, equation (1721) gives 
the well-known equation for elasticity of cylinder and pipe. 


If u = u r {r) e r is the displacement vector for non-fractal materials in 3-dimensional case 
(D = 3), then the strain tensor £p (r) has the following nonzero components that can be defined 
by 



E r r (G r , gr8.d ’ 

(73) 


£ w = div u - (e r , grad-u r ) = 

(74) 

and the invariant 

,. du r u r 

e = £ k k = div u = + . 

dr r 

(75) 


These equations with divergence and gradient can be generalized for non-integer dimensional 
case for u = u r (r ) e r . 

For non-integer dimensional model of fractal materials, we can use the definitions for the 
components of displacement vector in non-integer dimensional space in the form 


where we use the invariant 


As a result, we get 


£ rr — (e r , Grad r w r ), 
e w = Divf u - (e r , Gradf u r ), 


^ D du r D — 2 

e = £ kk = Div u = —— -| - u r 

dr r 

du r 

^rr ~X ) 

dr 


D -2 




-U r - 


(76) 

(77) 

(78) 

(79) 

(80) 
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Let us assume that the elastic constitutive relations for fractal material in isotropic case has 
the usual form 

U ij \skkSij T 2p J Sj t j. 

In this case, the non-zero components of stress tensor are 

du r D — 2 

u rr = 2 u £ rr + A e — (2a + A) — 7 -— T A-u r , 

or r 


<T, 


W 


_ . . du r D — 2 

— 2 /i + A e — A——|- (2/i + A) — 


-u r 


<9w r 


D -2 


^zz 2 /i £ 22 -(- A e A + A M r , 

ar r 


(81) 

(82) 

(83) 

(84) 


where we use £ 22 = 0. For D = 3, we have the usual constitutive relations for isotropic case in 
cylindrical coordinates. 


5.2 Elasticity of cylindrical fractal pipe with inside and outside pres¬ 
sures 

Let us consider the deformation of a fractal solid cylindrical pipe with internal radius R\ 
external radius R 2 with a inside pressure p\ and outside pressure P 2 . We use the cylindrical 
coordinates with the 2 -axis along the axis of the pipe. When the pressure is uniform along 
the pipe, the deformation is a purely radial displacement u = u r (r) e r , where e r = r/r. The 
equation for the displacement u r (r) in fractal pine is 


d 2 u r (r ) D — 2 du r 
dr 2 r dr 


D — 2 
—^— u, 


0, 


where 0 < D < 3. If D = 3, we get the usual (non-fractal) case. 

The general solution of equation (j85]h where D 7 ^ 1 and D 7 ^ 2, has the form 


(85) 


u r [r) = C\ r + C 2 r 2 ~ D . ( 86 ) 

Equations (j85|) with D — 1 has the general solution 

u r (r) = Ci r + C 2 r ln(r). (87) 

For D = 2, equations (1851) has the solution 

u r (r) = Ci + C 2 r. ( 88 ) 


Note that that fractal dimension of the pipe material can be D = 1 or D = 2. These cases do not 
correspond to the distribution of matter along the line and surface. These fractal dimensions 
describe a distribution of matter in 3-dimensional space (in the volume of pipe) such that the 
mass dimension is equal to D. 
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The constants C\ and (7 2 are determined by boundary conditions. Using that inside pressure 
is pi and outside pressure is P 2 , we get the boundary condition in the form 


Using (1861) . we get 


(J rr (R,\ ) Pi 5 (7 rr (R2) P 2 - 

du r 


dr 


— Ci + (2 — D ) C 2 r 


1 -D 


D -2 


u r = (D — 2) Ci + (D — 2) C 2 r 


l—D 


]~) — 2 

u rr = (2 p, + A) —— + A- u r = (2 fi + A (D — 1)) Ci — 2 p (D — 2) Co r 1 ~' D . 

or r 


Then 


The boundary condition (189|) has the form 

(2/i + A {D - 1)) Ci - 2 pi (D - 2) C 2 = -pi, 

(2/i + A (T> - 1)) Ci - 2 fi (D - 2) C 2 R\~ d = - p 2 . 

As a result, we have 

Ci = 


PiR\ D - p 2 Aj 


l-D 


>1 -D to1-D\ ’ 


(2/i + A (£> — 1)) (i^"" - 




P 2 “Pi 


jl-D nl-D\' 


2p(D-2) {R 2 -R 1 


The stress is 


Pi R\ ° ~ P2 R\ 




1 -D 


P2 ~P 1 


,1-D 


(^- D - f^) - /fi- 15 ) 

lf2<.D<3orl<.D<2, then we can rewrite equation (197|) in the form 




Pi Ri 1 - P 2 R 2 1 

(R?- 1 - R?- 1 ) 


P 2 ~ Pi 


(R 2- 1 - R?' 1 ) 


Ri Ro 


d- 1 


(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 


For the boundary conditions cr rr .(f? 2 ) = 0 and a rr (Ri) = —p, i.e. p 2 = 0 and pi = p for 
we have the solution of the form 




pRi 


D -1 




1 - 




D-P 


(99) 


This is the deformation of cylindrical pipe with a pressure p inside and no pressure outside. 
For (D = 3) we have 




pRi 


(Rl - RI) 


1-1 A 


( 100 ) 


that describes the stress of non-fractal material of pipe. 
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5.3 Rotating cylindrical fractal pipe 

Let us consider the deformation of a fractal solid that is describes by equation with external 
force f(r) for the displacement field u r (r) in fractal pine 


d 2 u r (r ) D — 2 du r D — 2 1 

—a"~2-' o- 2 — Ur + \ , o f ( 

or r or r z A + 2 ji 

where D > 0. The general solution of equation (1 101 1) has the form 


( 101 ) 


«r(r) = Cl r + C 2 r 2 D - (D - l)(X + 2 /M rdr ~ ° J f(r)r D 'dr^. 

Equations (IIOIH with D = 1 has the general solution 

r ( f Ri f R2 \ 

u r (r) — Ci r + C 2 r ln(r) + --( / f[r) ln(r) dr + ln(r) / f(r)dr\. 

A + 2 /r VJ Rl J Rl J 

For D = 2, equations (j 101 j) has the solution 


u r (r) = Ci +C 2 r 


1 

A + 2 n 




/(r)dr). 


( 102 ) 


(103) 


(104) 


Let us consider the deformation of a fractal solid cylindrical pipe with internal radius R± 
external radius R 2 rotating uniformly about its axis with angular velocity oj. Then the density 
of the centrifugal force is 

f r {r) = p 0 uj 2 r. (105) 

We use the cylindrical coordinates with the z-axis along the axis of the cylinder. When the 
pressure is uniform along the pipe, the deformation is a purely radial displacement u = u r {r ) e r , 
where e r = r/r. The equation for the displacement u r {r ) in fractal material is 


d 2 u r (r) D — 2 du r D — 2 p 0 u 2 

dr 2 r dr r 2 Ur A + 2p 

The general solution of equation (11061) has the form 


(106) 


u r {r ) = C\ r + C 2 r 2 D — A r 3 , 


(107) 


where 


A 


Po^ 2 

2(D + 1)(A + 2/i) 


(108) 


Using the condition that external forces do not act inside and outside the fractal pipe, we 
have the boundary condition 


a rr (Ri) = 0 , <J r r(.R 2 ) = 0 . 


(109) 
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Using (110711 . we get 

^ = Ci + (2 - D) C 2 r l ~ D - 3 A r 2 , 

Ur = (D-2)C 1 + (D- 2) C 2 r 1_D - A (D - 2) r 2 . 
r 

Then 

. . du r D — 2 

oy r = (2/i + A)—-h A-= 

or r 

= (2/i + A (£> - 1)) Ci — 2 (D — 2) C 2 r 1_D - 4 (6 /i + A (£> + 1)) r 2 . 
The boundary condition (I109H has the form 

(2/i + A (£> — 1)) Ci - 2 n {D - 2) Us R\~ d = A (6/i + A (D + 1)) 


( 110 ) 

( 111 ) 


( 112 ) 


(113) 


(2/i + A (D — 1)) Ci — 2 /i (£> - 2) C 2 i?2 _D = ^ (6 /1 + A (D + 1)) R 2 2 . 

Then 

A (6 M + A (D + 1)) (R? +1 - R? +1 ) 

1 (2/i + A( J D-l))(i? 2 D - 1 -i?f- 1 ) ’ 

r _ A (6 /i + A (D + 1)) (i? 2 - fl 2 ) (Ex l^- 1 
2 2fi(D-2)(R?~ 1 -R?- 1 ) 

Substitution of (111511 and (111611 into (110711 

^(6^ + A(£> + l))(«? +1 -< +1 )_, 

1 (2/i + X(D — 1)) (R 2 ~ l — -R?” 1 ) 


A(6n + X{D+1)) {Rl - R\) (R 1 R 2 ) D ~ l 
2 fi(D — 2) (R 2 ~ l — Ri _1 ) 

where A is defined by (I108|l . 

For the fractal cylinder (Ri = 0, R 2 = R), we have 


dr 3 . 


u r (r) 


Po u 2 ( 6/1 + A (D + 1) ^2 r _ r 3 

2(D + 1)(A + 2/i) y 2/i + A (id — 1) 


For D = 3, equation (II 18(1 gives 


u r (r) 


( 3 ^ + 2A R 2 _ 8 

8(A + 2fi) V A+ a 


(114) 

(115) 

(116) 


(117) 


(118) 


(119) 


Equation (111911 describes the displacement field for elastic cylinder with non-fractal material. 
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6 Gradient elasticity model for fractal materials 

6.1 Gradient elasticity theory of non-fractal material 

In the papers [63] EU [65] suggested to generalize the constitutive relations (1331) by the gradient 
modification that contains the Laplacian in the form 

(Tij T 2 fiEij'j lg A ^A £kk$ij ~t~ 2 flEij'j , (120) 

where l s is the scale parameter [[66] . 

For homogenous and isotropic materials equation for (11201) has the form 

A ^1 ± l 2 s Aj graddivu + 2/x ^1 ± l 2 s Aj Au + f = p D 2 u, (121) 

where f is the vector field of external force density. 

Using relation 

grad div u = curl curl u + Au, (122) 

we can rewrite equation (I12ip as 

A ^1 ± l 2 s A j curl curl u + (A + 2/x) ^1 ± l 2 s A j Au + f = p D 2 u. (123) 

If we assume that the displacement vector u is everywhere radial and it is a function of r = |r| 
alone (u^ = Wfc(|r|)), then 

curl u = 0. 

As a result, equation (1123|) has the form 

(A + 2p.) (l ±l 2 s A) Au + f = pD 2 u. (124) 

This is gradient elasticity equation for homogenous and isotropic materials with the spherical 
symmetry. For the non-gradient model (l 2 = 0), equation (I124|) for the displacement gives (TH61) . 

6.2 Gradient elasticity of fractal material 

A formal generalization of equations (11241) for fractal material in the framework of continuum 
models with non-iuteger dimensional space, where the displacement vector u = u(r, t), does 
not depend on the angles, has the form 

(A + 2 y) (l ± l 2 s (D ) V A^ V A? u + f = p D 2 u. (125) 

This is fractional gradient elasticity equation for homogenous and isotropic materials with the 
spherical symmetry. Let us consider equation (1125|) for static case (D 2 u = 0) with a minus in 
front of Laplacian, i.e. the GRADELA model for fractal materials [65] . 

(A + 2/i) (l - l 2 s (D) V A^ V A? u + f = 0. (126) 
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We can rewrite this equation as 


y Af ) 2 u - l~ 2 (D ) y Af u - (A + 2/r)" 1 l~ 2 (D ) f = 0. 


Using the vector Laplacian (HI) , we have 


d 2 u r (r ) D — 1 du r {r) D — 1 


y A?ufr) = ( ' 2 7 + 


n r (r)j e r . 


and 


y A D ) 2 u(r) = t^VCr) , 2(D-1) <9V(r) 


+ 


(£>-!)(£>-5) <9 2 i/ r (r) 3(U> - 1)(U> -3) du r (r) 3(D - 1)(D - 3) 


+ 


^>2 2 i ^3 

Using (11281 - fl29|) and f(r) = f(r) e r , equation (I127|) gives 

d 4 u r (r) 2(D - 1) d 3 u r (r) f (D — 1)(D — 5) 


n r (r)) e r . 


dr 4 


+ 


c>r 3 




<9 2 w r 


dr 2 


3(iA — 1)(D — 3) | ; _ 2 m D — 1^ du r {r) 


+ K ( D Y- 


dr 


3(0 3> + 17 2 (D)^ ) « r (r) — (A + 2/i) -1 /; 2 (D) /(r) = 0. 

General solution for the case f(r) = 0 is 


u r (r) = C 1 r + C 2 r l D - C 3 //(£>, r) - C 4 I K (D , r), 
where Ij(D,r ) and Ix{D,r) are the integrals of the Bessel functions 

Ii(D,r) — Dr j drr^ 0 ^ 1 j drr D ^ 2+1 I D / 2 (r/l s (D)), 

I k (D, r) = Dr J drr~ D - 1 j drr D / 2+1 K D/2 {r/l s (D)), 
where I a (x) and K a {x) are Bessel functions of the first and second kinds. 


(127) 

(128) 

(129) 


(130) 

(131) 

(132) 

(133) 


7 Thermoelasticity of fractal material 

Let us consider a generalization of thermoelasticity mmm for fractal material. In this sec¬ 
tion we consider a non-integer-dimensional model of thermoelasticity of fractal material. Note 
that thermoelasticity of fractal materials in the framework of the fractional-integral continuum 
model [EII91121], has been considered by Ostoja-Starzewski in Haim 
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7.1 Thermoelastic constitutive relation for fractal material 

If the isotropic material is non-uniformly heated, then the constitutive relation for a thermoe¬ 
lastic material must include [62] the term 

crij = A E k k Sij + 2 [i Eij - K a (T -T 0 ) S ij: (134) 

where A and /i are the Lame coefficients, K is the bulk modulus or modulus of compression. The 
third term in equation ()134p gives the additional stresses caused by the change in temperature. 
The bulk modulus for fractal materials is related to the Lame coefficients by 

K — A + —/i. (135) 

In this formula, we use the dimension D instead of 3 because 8kk = D for non-integer dimen¬ 
sional space (see Property 4 in Section 4.3 of [32]). 

If external forces being absent, then the stress is equal to zero c7{j = 0 and we have a free 
thermal expansion. Using a VJ = 0, equation 11134)1 gives 

A EkkSij + 2 iiEij - K a(T — T 0 ) S iy = 0. (136) 

Using 5kk = D and (1135)1 . we obtain 

£kk = Oi{T — T 0 ). (137) 

Because the function e = Ekk describes the relative change of volume caused by deformation, 
then a is the thermal expansion coefficient of the material [ 62] , 

The constitutive relation (1134)1 for fractal material can be represented in the form 

&ij = Ae kkSij + 2 iiEij - (D A + 2 /i) a T (T - T 0 ) 8^, (138) 

where we use the dimension D instead of 3 since K = (D A + 2 id)/D. 

7.2 Thermoelastic equations for fractal material 

For homogenous and isotropic non-fractal materials, we have the thermoelasticity equation 


A grad div u + 2/rAu + f — (3A + 2/i) a grad T = p D\ u, (139) 

where f is the external force density vector field. For the case u = u(r, t) and f = f(r, t), 
equation ()139)1 gives 

(A + 2/i) A u + f — (3A + 2/i) a grad T = p D\ u. (140) 

The thermoelasticity equation for fractal materials in the framework of models with non¬ 
integer dimensional space has the form 

(A + 2/i) v Af u(r, t) + f(r, t) — (D A + 2/i) a Gradf T(r, t) = pH t 2 u(r, t), (141) 
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where we assume f 
in the form 


f r (r, t ) e r and u = u r (r, t ) e, r . Using (J4J , equation (1 14 ID can be represented 


d 2 u r (r,t) D — 1 du r (r,t ) 

dr 2 r dr 


D - 1 


u r (r,t)+ 


-fr(r,t ) 


a(D A + 2/i) <9T(r, f) 


P 


D 2 u r (r , £). 


(142) 


A + 2/t'" v ' A + 2/i dr A + 2/a 

If the fractal material is non-uniformly heated, then the equation of equilibrium has the form 


<9V(r) D — 1 du r {r) _ D - 1 (DX + 2/i)a dT(r) 

dr 2 r dr r 2 Ur A + 2 /i <9r 


(143) 


This is thermoelasticity equation in spherical coordinates for pure radial deformation of fractal 
materials with fractal dimension D. For D — 3, we get the usual equation for thermoelasticity 
of solid ball [62]. 

General solution of equation (I143jl has the form 


<ar) = C 1 r + ^ + <§A±M£ (rT ( r) 


D 


r>D — 


Y / r rj Tirj dr ), 


(144) 


where Gi and C 2 are defined by boundary conditions. 


8 Conclusion 

In this paper, continuum models with non-integer dimensional space are suggested to describe 
isotropic fractal materials. A generalization of differential vector operators for non-integer 
dimensional space is proposed to describe elasticity of fractal materials in the framework of 
continuum models. The differential operators of first and second orders for non-integer dimen¬ 
sional space are suggested for rotatioually covariant scalar and vector functions. We consider 
applications for elasticity theory in the case of spherical and axial symmetries of the fractal 
material. Elastic properties of fractal hollow ball and cylindrical fractal pipe with inside and 
outside pressures, rotating cylindrical fractal clastic pipe are described. Equations for thermoe¬ 
lasticity and gradient elasticity of fractal materials are solved. 

Although the non-integer dimension does not reflect all properties of the fractal material, 
the suggested models with non-integer dimensional space nevertheless allows us to derive a 
number of important conclusions about the behavior of fractal materials. Therefore continuum 
models with non-integer dimensional spaces can be successfully used to describe elasticity and 
thermoelasticity of fractal materials. 

The proposed continuum models of fractal materials can be extended to more complex 
fractal materials: (1) We assume that continuum models with non-integer dimensional space 
can be generalized for anisotropic fractal materials [73]; (2) The non-integer dimensional models 
of fractal elastic materials can easily be generalized for the boundary dimensions d 7 ^ D — 1 
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[75]. (3) We also assume that differential and integral operators of fractional orders can also be 
defined for non-integer dimensional spaces to take into account non-locality of fractal materials. 
Note that dimensional continuation of the Riesz fractional integrals and derivatives [72] to 
generalize differential and integrals of fractional orders for non-integer dimensional space has 
been considered in [73]. 
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